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Abstract. Anomalous X-ray Pulsars (AXPs), thought to be magnetars, exhibit poorly understood deviations from a simple 
spin-down called "timing noise". AXP timing noise has strong low-frequency components which pose significant challenges 
for quantification. We describe a procedure for extracting two quantities of interest, the intensity and power spectral index of 
timing noise. We apply this procedure to timing data from three sources: a monitoring campaign of five AXPs, observations 
of five young pulsars, and the stable rotator PSR B1937+21. 
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INTRODUCTION 

Most pulsars spin with extreme regularity. Some have 
spin-down stability comparable to the best atomic 
clocks. Others, however, appear to deviate from a simple 
spin-down, requiring a polynomial of high degree to 
model their rotation over even a year (see Figure Q] 
for some examples). Still others are sufficiently unpre- 
dictable that phase-connecting the observations — that 
is, exactly counting the number of revolutions between 
observations — is difficult or impossible. The causes of 
such deviations, called "timing noise", remain poorly 
understood IjJ] . 

Anomalous X-ray pulsars (AXPs) are slow pulsars 
(periods 2-12 s) whose X-ray luminosity exceeds the 
power available from spin-down. The leading description 
of these objects is the magnetar model [2]. Eight AXPs 
are known, plus two candidates^ AXPs exhibit a tremen- 
dous range of variability in their radiative behaviour, in- 
cluding flux and profile variations (01, |01), bursts JH], 
glitches [6], and strong timing noise. In fact, the AXP 
IE 1048.1—5937 has intervals during which phase con- 
nection is difficult to impossible and periods where it is 
straightforward (in spite of an unchanged cadence of ob- 
servations). This indicates that the nature or amplitude 
of its timing noise changes. By quantifying the timing 
noise, we hope not only to describe these changes but to 
probe the differences between AXPs and other classes of 
pulsar. 



1 An online catalog of AXPs and SGRs is maintained by Cindy Tarn at 

http : / /www . physics . mcgill . ca/~pulsar /magnetar /main 
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FIGURE 1. Timing noise. Observations of the AXP IE 
1048. 1 —5937 (top left) are compared with algorithmically gen- 
erated timing noise with power spectral indices of —4 (top 
right), —2 (bottom left) and (bottom right). In all cases the 
best-fit cubic has been subtracted (accounting for two deriva- 
tives of frequency). Observe that timing noise with a large neg- 
ative spectral index takes on a smooth appearance because of 
the strong low-frequency components. 



ANALYSIS 

The analysis of timing noise poses a number of techni- 
cal difficulties. Observed timing noise has a steep red 
spectrum, with power spectral indices often near —4 
(see Table [TJ. Noise with such a red spectrum cannot 
be analyzed with normal Fourier techniques, as individ- 
ual Fourier bins do not reject low frequencies strongly 
enough (this is often called "spectral bleeding") |7j]. Fur- 
ijhgy/iie span of data limits spectral resolution at low fre- 
quencies no matter what technique is used. 

To bypass the limitations of Fourier analysis, we are 
developing a method that we believe will improve spec- 
tral analysis of red noise. Our code constructs a filter- 
bank for each data set based on orthogonal polynomials, 
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FIGURE 2. Spectral responses of orthogonal polynomials 
and periodogram bins, sampled at the times of our observations 
of IE 1048.1—5937. The top plot compares the response of a 
degree 9 orthogonal polynomial and a periodogram bin to a 
sine wave. The bottom plot shows the same response scaled 
by /~ 4 . This indicates which parts of a red spectrum the 
niters respond to. Observe that the periodogram is much more 
sensitive to the lowest frequencies in such a red spectrum than 
it is to the frequencies at its nominal peak, while the orthogonal 
polynomial rejects them strongly. 
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FIGURE 3. A model fit to observations of IE 1048.1-5937. 
The circles are the amplitudes of the responses of an orthogonal 
filterbank to the residuals pictured in Figure [T] The horizontal 
error bars indicate the full width at half-maximum response of 
each filter. If the best-fit model we computed is the true model, 
we expect the amplitude returned by each filter to be distributed 
according to a X 2 distribution with one degree of freedom. The 
line connects the expected value for each such distribution, and 
the vertical error bars attached to it indicate the ±1(7 points. 
The horizontal line indicates the background level of white 
noise (coming mostly from observational uncertainties); this is 
also indicated with vertical error bars attached to the circles. 
For most of the plotted filters the observational uncertainty 
is negligible compared to the variability we expect from the 
timing noise. 



extending the procedure described by Deeter and Boyn- 
ton [0. Specifically, Deeter and Boynton [7] describe a 
procedure for computing a family of polynomials orthog- 
onal with respect to the white observational noise in a 
data set. The resulting filters reject low frequencies very 
strongly (see Figure |2). They then carry out a second or- 
thogonalization procedure based on a particular red noise 
model, and they select a single filter from the resulting 
filterbank.To examine higher frequencies, they choose 
to apply this process from scratch on halves, quarters, 
eighths, etc. of the original time series. 

By comparison, we compute a family of polynomials 
orthogonal with respect to the white noise, as Deeter and 
Boynton |7[] do. However, we simply use these polyno- 
mials as our filterbank; their frequency responses vary, 
but as a filterbank they provide reasonable frequency res- 
olution at high frequencies, and strong rejection of low 
frequencies. The polynomials remain covariant with re- 
spect to timing noise, so we take that into account in 
our fitting procedure. Given a power-law model of timing 
noise (as described in Scott et al. ||8|]), we can compute 
the likelihood of the observed filter values using the fil- 
terbank's covariance matrix with respect to this red noise 
model. We then carry out a computationally expensive 
numerical optimization procedure to find the maximum- 
likelihood power-law model. 

The power-law spectral model we fit to the timing 
noise is 




where F is the power spectral density in s 2 /Hz and to is 
the frequency of the timing noise. Fitting then gives us 
A and a. We also take into account the uncertainties on 
our arrival-time measurements, which are independent. 
This white timing noise, due to observational effects, sets 
an upper limit on the frequencies we can include in our 
model. 



DATA 

We are considering monitoring data for three kinds of 
sources: five AXPs, five young pulsars, and the stable 
rotator PSRB 1937+21. 

Of the known AXPs, we consider monitoring data 
five (IE 1048.1-5937, IE 2259+586, 4U 0142+61, IE 
1841-045, and RXS J170849.0-400910) on a regular 
basis using the Rossi X-ray Timing Explorer (RXTE), 
as described in Dib et al. 2007 (this volume). Observa- 
tions, spanning a decade, were generally carried out on a 
weekly or monthly basis, and yielded time-of-arrival un- 
certainties generally less than 1 % of a period. These data 
sets have been phase-connected where possible. Glitches 
have been observed in all these AXPs. We consider only 
inter-glitch intervals, as it is frequently impossible to un- 
ambiguously count phases across a glitch, and because 
glitches seem likely to be produced by a different mech- 
anism than timing noise H. We were also forced to ex- 
clude intervals where the pulse profile changed signifi- 



TABLE 1. Preliminary timing noise model fits. Each row represents a single span of data; separate 
spans for the same source were fit independently. 



Source 


Exponent 


log amplitude (s 2 /Hz) 


Period (s) 


Dates (MJD) 


IE 1048.1-5937 


-4.3±0.2 


10.8±0.3 


6.5 


53030-54111 


IE 1841-045 


-4.3 ±0.9 


8.9±0.4 


11.8 


51225-52438 




4.3 ±1.1 


9.1±0.7 




53030-53816 


1RXS J170849.0-400910 


-5.7±2.3 


8.1 ±2.3 


11.0 


53638-54016 


4U 0142+61 


-4.4±1.3 


6.3±0.4 


8.7 


51611-53787 


IE 2259+586^ 


-2±2 


4.7 ±0.9 


7.0 


53012-54181 


PSR B0540-69 


-4.7±0.8 


3.2±0.4 


0.050 


51342-52936 


AX J1811.5— 2134 


-3.7±1.3 


3.4±0.6 


0.065 


52342-53333 


PSR B 1937+21 


-3.6±1.2 


-6.6±0.3 


0.0016 


45986-49303 



* The timing noise model is a poor fit to the data for this source 



cantly (making identification of a fiducial point difficult), 
intervals where phase-connecting the observations was 
impossible, and intervals where the span of good data 
was too short to be useful. 

We also consider a collection of observations (pri- 
marily using RXTE) of five young pulsars. In particular: 
we consider 7.6 years of observations (interrupted by a 
glitch) of PSR B0540-69 [ 10]; we consider 2.7 years of 
observations of PSR J 1 8 1 1 — 1 925 (in the supernova rem- 
nant (SNR) Gl 1.2— 0.3); we consider 7 years of observa- 
tions (interrupted by a data gap) of PSR J1846-0258 (in 
the SNR Kes 75) 111 ill ; we consider 2.7 years of combined 
X-ray and radio (using the Robert C. Byrd telescope at 
the Green Bank observatory and the Lovell telescope at 
Jodrell Bank) observations of PSR J0205+6449 (in the 
SNR 3C 58) (Livingstone et al. 2007, these proceedings); 
and we consider 21.3 years of combined X-ray and radio 
observations (with the Molonglo Observatory Synthesis 
Telescope and Parkes) of PSR B1509-58 lUl . 

The pulsar PSR B 1937+2 1 is a millisecond pulsar with 
exceptionally stable rotation. For this reason its rotation 
has been monitored carefully (see, for example, Kaspi 
et al. II 1310 . and we consider 9.0 years of radio observa- 
tions obtained with the radio telescope at Arecibo. 

RESULTS AND DISCUSSION 

Our code still requires work. Nevertheless we can give 
some preliminary measurements that make quantitative 
the claim that AXPs have very strong timing noise; see 
Table Q] Observe that the AXP timing noise is, for the 
most part, vastly stronger than the timing noise in the 
other sources, and although all the exponents are consis- 
tent with —4 (the value for white torque noise) there is 
a hint that AXPs may have slightly redder timing noise 
than other pulsars. 

We hope to be able to use our code and data to answer 
several questions. Do AXPs have time- varying timing 



noise? Does the timing noise in AXPs differ in character 
from the timing noise in other classes of pulsar? Can we 
estimate the degree to which measurements of proper 
motions, braking indices, and so on are contaminated 
by timing noise? We also hope eventually to extend our 
analysis to a broader range of radio pulsars, and possibly 
also to Soft Gamma Repeaters, if this is feasible. 
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